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ABSTRACT 

This paper presents the full VLT/VIMOS-IFU data set and related data products from an 
ESO Large Programme with the observational goal of obtaining two-dimensional kinematic 
data of early-type lens galaxies, out to one effective radius. The sample consists of 17 
early-type galaxies (ETG) selected from the SLACS gravitational-lens survey. The galaxies 
cover the redshift range from 0.08 to 0.35 and have stellar velocity dispersions between 200 
and 350 km/s. This programme is complemented by a similar observational programme on 
Keck, using long-slit spectroscopy. In combination with multi-band imaging data, the kine- 
matic data provide stringent constraints on the inner mass profiles of ETGs beyond the lo- 
cal universe. Our Large Programme thus extends studies of nearby early-type galaxies (e.g. 
SAURON/ ATLAS 3D) by an order of magnitude in distance and toward higher masses. We 
provide an overview of our observational strategy, the data products (luminosity-weighted 
spectra and Hubble Space Telescope images) and derived products (i.e. two-dimensional fields 
of velocity dispersions and streaming motions) that have been used in a number of published 
and forthcoming lensing, kinematic and stellar-population studies. These studies also pave 
the way for future studies of early-type galaxies at z ~ 1 with the upcoming extremely large 
telescopes. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: structure — galaxies: kine- 
matics and dynamics — techniques: spectroscopic — gravitational lensing 



1 INTRODUCTION 

The formation and evolution of early-type galaxies (ETGs) is 
a widely studied topic in present-day astrophysics, in partic- 
ular due to a number of tight correlations between their ob- 
servables, such as structural parameters, kinematics, colours 
and c entral supermassive black-hol e masses (e.g . Illingworth 
1977; Sandage & Visvanathan 1978; Diorgovski & Davis 1983; 
Dress ier et al. 198 7; Ferrarese & Merritt 2000; Gebhardt et al] 
200Qh . SeelE enzm I (l2006h for a recent review. These correlations 



tion paradigm fe.g. 'Binnev"l978': 'Blum enthal et al.| ll984'). where 
ETGs are seen as the end products of the st ochastical merging 
of sn i aller bui lding blocks (e.g. spiral galaxies: [Toomre & Toomrd 
ll972l : lBarnesl[T99 2). 

Both observational and theoretical studies have progressed 
rapidly over the last decade, largely thanks to major observational 
programmes conducted on the ground and from space (e.g. 2dF, 



seem hard to explain at first sight in the hierarchical CDM forma- 
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programmes conauctea on tne grouna ana irom space (e.g. zar, 
SDSS COSMOS, COMB017;IColless etal]l200ll : lBernardi et all 
l2003l : IScarlata e t al.ll2007l : iFaber et al.ll200 7') and due to the large- 
scale numeric al simulations now adays affordable by supercom- 
puters (e.g. iHoDkins et al.l l2006h . The SAURON and ATLAS^^ 
projects provide detailed analysis on substantial samples of nearby 
ETGs based on deep integral-field spectr oscopic data and serve 
as a local benchmark for ETG studies (e.g. ICappellari et al.ll2007l : 
Emsellem et al. 2007, 2011). 

Despite the tremendous increase in data volume and computa- 
tional power, many questions remain very difficult to answer. This 
is due to several factors: (i) the types of available data (usually only 
luminosity- weighted kinematics and imaging) and their quality (in 
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terms of signal-to-noise ratio and spatial/spectral resolution) have 
remained limited compared to their volume, (ii) mass modeling 
techniques often suffer from intrinsic degeneracies that are hard 
to overcome without very high quality data that can only be ob- 
tained for small samples of nearby ETGs. ETG studies that are now 
conducted beyond the local Universe with 8-lOm-class and space- 
based telescopes suffer from many of the same issues that similar 
studies of the local Universe, with 4 m-class telescopes, faced more 
than a decade ago. 



One of the prevailing and, even today, still not precisely 
answered questions is how important dark n iatter is in the in- 



ner barvon-dominated redons of ETGs (e.g. ISagliaetal. 


1992|, 


1993': 'Bertola et al." 


1993"; "Bertin et al.""l994; Carollo et al. 


1995|; 


Gerhard et al.. 1998 


: Loewenstein & White .1999: Gerhan 


i et al. 


20011: iKeetonI 12001 


: Romanowskv et al.l 20031: iMamon & 


:.okas 


20051: Ide Lorenzi et al.l 120091). Open questions include in particu- 



lar the following: (i) How much dark matter is there precisely in- 
side the inner few effective radii of ETGs and how is it distributed? 
(ii) How do the stellar and dark matter distributions of ETGs evolve 
with cosmic time and are there trends with galaxy mass? (iii) Do 
these observations agree with theoretical predictions? 



Because dark matter is not directly observable but can only 
be inferred from other observations, these questions are particu- 
larly difficult to answer for more distant ETGs where data qual- 
ity progressively deteriorates due to lack of signal-to-noise, even 
with the largest ground-based telescopes. To overcome a num- 
ber of these hurdles, several systematic programmes were initi- 
ated over the last decade to combine the constraints of strong 
gravitational lensing with those of stellar kinematics (LSD and 
SLAGS; e.g . Koopmans & Treu 2002; Treu & Koopmans 2004; 
Bolton et al. [20061: ITreu et al. 2006; Koop mans et aLlOOOd 120091 : 
Bolton et al.l2008l : lAuger et al. 2009). This combination has turned 
out particularly powerful in breaking degeneracies in ETG mass 
models, despite the limited quality of the kinematic da ta (e.g. 
Treu & Koopmans 2004; Koopmans et al. 2006, 2009; Auge r et aP 
2010). The reason is that the mass enclosed by the Einstein radius 
of a lens galaxy can be accurately determined and breaks the usual 
mass-anisotropy degeneracy of ETGs to a large extent (see e.g. 
iKoopmansI |2004 for an explanation based on a simple power-law 
model). 

Whereas the LSD and SLAGS programmes used predomi- 
nantly luminosity-weighted stellar velocity di spersions and were 
modeled using the spherical Je ans equation (e.g. lKoopmans & Treul 
l2003l : [ireu & Koopmans"2004), one might argue that (slowly) ro- 
tating systems and/or systems with orbital anisotropy cannot be 
modeled prop erly that way (see for example the discussion in 
lKochanel j'2006). To address these valid issues, a new and substan- 
tially larger observational programme of a sub- sample of SLAGS 
lenses, using obse rvations with the VLT and Kec k, was started in 
2006 (.Czoske et al.l2008lJ2009l : rBarnabe et al.l2010) to obtain two- 
dimensional kinematic fields (both first and second moments of the 
line-of- sight velocity distribution) out to their effective radii, com- 
plemented with a rigor ous modeling effort based on two-integral 
axisymmetric models ( Barnabe & Koopmansll2007l : lBarnabe et al.l 
12009) . The self-consistent combination of gravitational-lensing and 
stellar-kinematic data sets allows better constraints to be set on the 
mass distribution of ETGs ( Czoske et a l.l2008l : lBarnabe et al.l2009l . 
llOlO, 2011), but also to test less rigorous methods that might still 
have to be used at even higher redshifts (i.e. z > 0.5) until IFU ob- 
servations from the next generation of telescopes (such as ESO's 
ELT) become available. 

Whereas the above-mentioned results show that these data 



provide more precise and accurate results on the total mass pro- 
file of ETGs in the inner few effective radii, several problems still 
remain that are critical for our understanding of galaxy formation. 
For example, what fraction of the mass inside several effective 
radii is due to dark matter and what are the relative mass distri- 
butions of the stellar and dark-matter components? This is much 
harder to answer with kinematic and lensing data alone, since there 
are many combinations of the stellar and dark-matter distributions 
that lead to similar lensing and kinematic data (Treu & KoopmansI 
l2004l : lBarnabe et allllOOT) . If, on the other hand, the stellar mass- 
to-light ratio can be determined independently, these degeneracies 
can mostly be lifted. This assumes, however, that the broad-band 
colours provide sufficient information to break degeneracies in the 
IMF models, which in general they do not. Still, current constraints 
can already set limits on the IM F and stellar mas s -to-light ra- 
tios that are becoming competitive (iGriUo et al.ll2008l : lAuger etaP 
2009, 2010; Spiniello et al. 201 ll). To go beyond these broad-band 
inferred stellar mass-to-light ratios, one can also use the same (IFU) 
spectra to determine line indices from which more precise ages, 
metallicities and mas s-to-light ratios of the ETG stellar population 
can be obtained fe.g. lTrager et al.|[200 Q^. Whereas this information 
is going to be used in forthcoming publications, here we concen- 
trate on describing the data for kinematic purposes. 

In this paper, we present the full VLTA^IMOS-IFU data of 
our sample of 17 SLAGS early- type lens galaxies, collected for 
these purposes. The sample is presented in Sect.|2] We describe the 
integral-field spectroscopic observations using VIMOS on the VLT 
in Sect. l3.1l and summarize the HST observations as needed for this 
project in Sect. 13.21 The data reduction of the IFU data is described 
in some detail in Sect. ID and the kinematic analysis, i.e. the deriva- 
tion of two-dimensional maps of line-of- sight velocity and velocity 
dispersion, in Sect.|5] In Sect.[6l we provide a list of additional ob- 
jects that are visible in the VIMOS/IFU field of view around the 
lens galaxies. Finally, we conclude in Sect. |7] Appendix lAl gives 
a recipe for estimating the expected photon noise on the reduced 
spectra based on a simple model. 

The fully self-consistent combined lensing/dynamics analysis 
of the d ata described in this p aper is presented in the companion 
paper bv lBarnabe et al ] (l201lh . 



2 SAMPLE 

The sample described in this paper consists of seventeen early-type 
galaxies from the Sloan Digital Sky Survey that have been con- 
firmed as secure gravitational-lens systems in the SLAGS surveyQ 
The SLAGS survey searches for gravitational lens systems in 
the spectroscopic data of a subsample of SDSS galaxies comprising 
the luminous red galaxy sample dEisenstein et al.ll20()ll) and pas- 
sive members of the MAIN sample lstrauss et al.ll2002l) . defined as 
having rest-frame Ka equivalent widths less than 1.5 A. The pres- 
ence of emission lines at a redshift higher than that of the target 
galaxy in the SDSS spectra indicates that there is a background 
galaxy within the fibre, which may have been lensed. Candidates 
were selected for follow-up snapshot imaging with HST, allow- 
ing in most cases a robust confirmatio n or rejection of the lens- 
ing hypothesis. See Bolton et al ] (12006) for details on the selection 
procedure. .Treu et al.. 1200^ conducted a number of tests to verify 



^ In fact, one system, J1250B, was downgraded to grade B ("possible lens") 
in lBolton et Z[ ( l2008h . 
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Figure 1. Distribution of the SLACS/IFU sample in redshift and velocity 
dispersion. The velocity dispersion measurements are taken from this paper, 
the error bars are stati stical errors. For comparison, the SLACS early-type 
sample of'Auser et al." ('2009, grey points), the SAURON sample at redshift 
z ^ (triangles: Emsellem et al. 2004) and a sample of Coma cluster early- 
type galaxies (stars: iMehlert et al.ll200Ql : iCorsini et al.ll2008D are shown. 

that the SLACS lens sample is statistically indistinguishable from 
the parent samples. Results on the structure of the early-type lens 
galaxies can therefore be taken as representative for the population 
of (massive) early- type galaxies as a whole. 

In order to be observable from Paranal, lens systems close to 
the equator (6 < 15°) were selected from the (mostly northern) 
SLACS catalogue. The final sample of 17 systems chosen for VI- 
MOS/IFU follow-up is evenly distributed in right ascension. 

The location of the 84 early-type SLACS galaxies studied by 
lAuger et all (l2009l) in the parameter space of redshift and velocity 
dispersion is shown in Fig. [T] The seventeen galaxies from the VI- 
MOS/IFU subsample are marked; they are representative of the full 
SLACS early-type galaxy sample. Fig. [T] also shows the 48 local 
early- type galaxies studied by the SAURON team (Emsellem et al.l 
I2OO4I), as well as 17 ear ly-type galaxies in the Coma cluster studied 
bv lThomas eTaP (l2007h . The VIMOS/IFU sample represents a ma- 
jor step out to cosmological distances compared to the SAURON 
and Coma samples. In terms of velocity dispersion, it overlaps with 
the more massive half of the SAURON sample and thus lends itself 
to comparison of the structure of early-type galaxies at z > 0. 1 to 
their local counterparts. 



3 OBSERVATIONS 

3.1 Integral-field spectroscopy 

We have obtained integral-field spectroscopic observations of 

17 systems using the integral-field unit (IFU) of VIMOS 

JEeFevre et al.ll200lh . 

VIMOS is a wide-field imager, multi-object and integral-field 
spectrograph mounted at the Nasmyth focus B of the Very Large 
Telescope UTS (Melipal) on Paranal, Chile. The integral-field unit 
samples the focal plane with a total of 6400 lenslets, each of which 
is coupled to an optical fibre. The opposite ends of the fibres are 
arranged to form pseudo- slits of 400 fibres each. The light emerg- 
ing from the pseudo- slit is dispersed in the usual way using grisms 
and optionally order- separating filters, and the resulting spectra are 



recorded on four 2k x 4k EEV CCD chips. In low-resolution mode, 
the spectra from four pseudo- slits are stacked in the dispersion di- 
rection on each of the four CCDs. The arrangement is such that 
each CCD records the spectra from one quadrant of the field of 
the IFU head. In high-resolution mode, only the spectra from one 
pseudo-slit fit onto a CCD and only the central 1600 elements of the 
IFU head are used. The spatial scale in the focal plane can be cho- 
sen to be 0.67 arcsec or 0.33 arcsec per spectral element. We used 
the large spatial scale and high-resolution mode and thus work with 
a field of view of 27" x 27" covered by 40 x 40 spectra. 

Three of the 17 systems of the sample were observed as ESO 
Programme 075. B -0226 (henceforth the "pilot programme") be- 
tween June 2005 and January 2006. The remaining fourteen sys- 
tems were the targets of an ESO Large Programme, 177.B-0682 
(the "main programme"), and were observed between April 2006 
and March 2007. For the pilot programme, the high-resolution blue 
grism (HR-Blue) was used, which has a nominal resolution of 
R = 2550 and covers the wavelength range 4000 A to 6200 A; for 
the main programme the HR-Orange grism was used, which covers 
a somewhat redder wavelength range between 5000 A and 7000 A 
at comparable resolution (R = 2650). 

All observations were carried out in service mode. For this, 
the total observing time was broken down in observing blocks 
(OBs), each with an execution time of one hour and comprising 
science and calibration exposures. For the pilot programme, three 
dithered science exposures of integration time 555 seconds on one 
target were taken during an OB. For the main programme, a single 
exposure per OB was taken, with integration time 2060 seconds. 
The main observational constraint for scheduling the OBs for this 
project was a maximum seeing of 0.8 arcsec FWHM (full width at 
half maximum). 

We followed ESO's standard calibration plan which comprises 
three flat-field exposures and one arc-lamp exposure per OB, taken 
immediately after the science exposures. Bias exposures were taken 
during day time in blocks of five exposures. 

Standard stars are observed at varying intervals by ESO staff" 
and were therefore in most cases not available for the same nights 
that our science data were taken. They are useful for relative flux 
calibration but not for absolute spectro-photometric calibration. 

3.2 HST imaging 

Modelling stellar dynamics and the lens efl'ect requires photomet- 
ric information on both the lens galaxy and the lensed background 
source. We use various HST observations obtaine d for the SLACS 
project (see lBolton et aLll2008l : lAuger et aLll2009h . 

Lens modelling aims at reconstructing the two-dimensional 
surface-brightness distribution of the background source. This is 
obtained from high-resolution ACS images from whic h elliptical B- 
spline models of the lens galaxies were subtracted f see lBolton et alj 
I2OO8I . for details). For eight systems, we have deep F814W images 
from Programmes 10494 and 10798 (Cycle 14, PI: L. Koopmans) 
consisting of four exposures. For the remaining systems only sin- 
gle exposure images from snapshot Programme 10174 (Cycles 13, 
PI: L. Koopmans) are available. The pixel scale of these images is 
0.05 arcsec. 

The surface-brightness distributions of the lens galaxies are 
obtained from NICMOS F160W images from Programmes 10494, 
10798 and 11202 (Cycle 16, PI: Koopmans). The full resolution 
images have a pixel scale of 0.05 arcsec and include four exposures 
for a total exposure time of 2560 s. In order to have proper flux 
weights for the kinematic maps we resample the NICMOS images 
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Table 1. The VIMOS/IFU sample. Apparent brightness my, effective radius /?efF v and Einstein radius %nst are from'Auser et al IbOOQ). Redshifts and SDSS 
velocity dispersion are from Bolton et al. (2008). VIMOS velocity dispersions were measured from aperture-integrated spectra from the IFU data; errors 
include only the effect of noise, not of template mismatch (5-10 %, Sect. 15. IV The last two columns specify the grism used for each object and the number of 
observing blocks spent on each target. 
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to the VIMOS/IFU grid with 0.67 arcsec per spatial pixel after con- 
volution to match the point spread function (PSF) of the VIMOS 
data, taken to be Gaussian with FWHM of 0.8 arcsec. 

Three systems (J 1251, J 1627 and J2321) were not observed 
with NICMOS. For these systems we use the B- spline models to 
the F814W ACS images to estimate the surface brightness of the 
lens galaxy. We add random noise according to the ACS weight 
maps to mimic observational data. 

All images were astrometrically registered using the measured 
position of the lens galaxy and in some cases additional objects 
visible in the VIMOS field of view. 

The F814W images of the systems in the sample are shown in 

Fig.ia 



4 DATA REDUCTION 

The VLT VIMOS/IFU data were reduced using the vipgi pack- 
age which was developed within the framework of the VIRMOS 
consortiu m and the VVDS project. vipgi H has been d escribed in 
detail bv IScodeggio et all (120051) and IZanichelli et alJ (|2005). In 
ICzoske etal! dlOOsL we have already given an overview of the data 
reduction procedure as applied to our data set, so that we restrict 
ourselves here to giving the main parameters. 

The wavelength calibration is based on about 20 helium 
and neon lines spread over the wavelength range. Fitting a third- 
order polynomial results in residuals with root-mean- square of 
~ 0.075 A, corresponding to ~ 5kms~\ significantly smaller than 
the spectral resolution and negligible in our analysis. 



http://cosmos.iasf-milano.inaf.it/pandora/vipgi.html 



After rectification of the two-dimensional spectra onto a lin- 
ear wavelength grid with a dispersion of 0.644 A per pixel, one- 
dimensional spectra a re extracted fo r each fiber in an optimal way 
using the algorithm o f lHornel(ll986h . We refer to the resulting col- 
lection of one-dimensional spectra as a "data cube" even though the 
data structure is, strictly speaking, not a three-dimensional cube. 
Spectra are labeled by the spatial coordinates L and M (or world 
coordinates a and 6). 

A relative flux calibration (giving correct flux as a function of 
wavelength) is done using the standard observation nearest in time 
to the execution of a given science observing block. Since this does 
not necessarily come from the same night as the science exposures 
and no attempt is made to observe them at the same airmass, the 
standard observations can only provide a relative correction but no 
absolute spectro-photometric calibration. 

The varying transmissivity of the fibres is corrected by mea- 
suring and comparing the flux under a sky line in the diff'erent ex- 
tracted spectra. We used [O i] A5577A for this purpose, integrating 
over a range of seven pixels (4.2 A) and subtracting an estimate of 
the underlying continuum from a window of five pixels on one side 
of the line. 

The observed galaxies are smaller than the field of view of 
VIMOS/IFU, leaving a sufficient number of fibres pointing at blank 
sky to permit good sky subtraction. These sky fibres are identified 
automatically from a histogram of the total fluxes recorded in each 
fibre, integrated over the entire wavelength range. The sky fibres 
are then grouped according to the shape of the [O i] A5511 line, a 
mean over the sky spectra in each group is computed and subtracted 
from each row according to the group it is judged to belong to. The 
sky subtraction works well even redwards of 6200 A, where the sky 
is dominated by the OH molecular bands. 

The data cubes from all the exposures on a field are finally 



Two-dimensional kinematics ofSLACS lenses -IV. 5 





J1251 




J1330 








f 





Figure 2. HST/F814W images of the systems in the sample. In images obtained in a snapshot programme cosmic ray hits were interpolated over for display 
purposes. 



combined into a final data cube by taking the median. Telescope 
off'sets between the exposures are corrected through spatial shift- 
ing by an integer number of fibres. Subpixel shifts could in prin- 
ciple be corrected for by interpolation between adjacent fibres but 
in practice it turns out that the centroids of the galaxies in image 



reconstructions of the individual exposures are hard to measure to 
subpixel accuracy, given the level of accuracy of the fibre transmis- 
sion correction and the resulting noise in the image reconstructions. 

A short-coming of vipgi is that it does not provide noise esti- 
mates on the reduced spectra. Since proper weighting is desirable 
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Table 2. Stellar templates used for kinematic analysis of the VIMOS/IFU 
sample 



Galaxy 


Template star 


Spectral type 


J 0037 


HD 249 


Kl IV 


J OZiO 




K2 III 


J0912 


n I J iz,i iH-vj 


K3 IV 


J 0935 


HD 195506 


K2III 


J 0959 


HD 102328 


K3 III 


J 1204 


HD 221 148 


K3 III 


J 1250A 


HD 145328 


Kim 


J 1250B 


HD 216640 


Kl III 


J 1251 


HD 145328 


Kim 


J 1330 
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for the determination of kinematic and stellar population param- 
eters, we reconstruct noise spectra from a simple model includ- 
ing photon noise and read-out noise. We require noise spectra for 
each spectral element in the data cube and for the aperture- summed 
spectra. It turns out that at each step of the reduction procedure, the 
noise estimate can be written as a rescaled version of the current 
intermediate data product, including the sky background. The nec- 
essary noise spectra can thus easily be produced by creating and 
rescaling a second data set, reduced in the normal manner but with 
sky subtraction turned off. A detailed derivation of the noise esti- 
mation procedure is given in Appendix [A] which can also be read 
as a walk-through of the data reduction process. 

Fig.[3]shows integrated spectra obtained by summing the one- 
dimensional spectra from fibres within elliptical apertures follow- 
ing the shape and size of each lens galaxy. Typical diameters of 
these apertures are around 4 arcsec, i.e. a little larger than the 3 arc- 
sec circular aperture used by the SDSS. 

The spectra are typical for early-type galaxies, with strong 
metal absorption lines, in particular Cai H and K, the G and the 
Mgb bands. With the exception of J2321 (Czoske et al. 2008), no 
significant emission lines from the lens galaxies are detected. This 
is true even for those galaxies that show clear disk and spiral struc- 
ture in the HST images. 

Most of the spectra clearly show [O ii] A3121 emission from 
the lensed source, although in some cases this line falls outside the 
wavelength range of the VIMOS spectra. In some cases, in par- 
ticular J1250B and J 1451, several B aimer emission lines from the 
source can be seen. These become more prominent in the residuals 
of the kinematic fits, described in Sect. (5] 

5 KINEMATIC ANALYSIS 
5.1 Method 

Many methods to determine the line-of- sight velocity distribu- 
tion (LOSVD) of earl y-type galaxies from spectra have been pro- 
posed in the past (e.g. |R ix & W hitdl 19921 : Ivan der Marel & Franxl 
ll993l : ICapDellari & Em sellem 2o6i). We use the conceptually sim- 
plest method, template-fitting in pixel space. Our method is im- 



plemented as a package (slacR) in the statistical language R 
(Ir Development Core Team 2007 fl following descri ptions from 
several auth ors, in particular Ivan der Marel & Franx (1993) and 
iKelson et"al . (2000). Compared to the version used in Czoske et al.l 
we have modified and extended our method sufficiently to 
warrant a full description of the method here. 

The model parameters 6 are determined by minimizing the 
merit function 

S =Yj^ist-s{A„0)f , (1) 

i i 

where the sum extends over all "good" pixels Ai in the observed 
spectrum Si. The observed spectrum is transformed to the rest frame 
of the respective lens galaxy prior to the analysis by dividing ob- 
served wavelengths by a factor 1 + ziens- The eff'ective resolution of 
the spectrum is increased to /?(1 + ziens) compared to the nominal 
resolution R of the spectrograph (^Sect. l3.ll l. 

The data-model s (in vector-matrix notation) is given by the 
equation 

Si = [G * t](Ai) p^'^\Ai) + q^''\Ai) + 6/ . (2) 

Here, the vector G is the LOSVD, taken to be a Gaussian in veloc- 
ity, with kinematic parameters v (streaming motion) and cr^ (veloc- 
ity dispersion). The vector ^ is a stellar template spectrum, taken 
from the Indo-US library of stellar spectra ( Valdes et al. 200i). 
Templates were chosen by fitting a random sample of IndoUS spec- 
tra to the aperture-integrated VIMOS/IFU spectra and selecting one 
of the best-fitting (in the least- squares sense) template candidates. 
We did not always choose the template giving the absolute mini- 
mum It was noticed that HD 195506 appeared among the best- 
fitting templates for several systems and was therefore chosen for 
all of these. Table |2] lists the templates used for the kinematic fits, 
as well as their spectral type and luminosity class. As expected, 
the selected template stars are predominantly late- type giants. The 
same template was used for fitting all the individual fibre spectra 
for each system. The eff'ect of template mismatch on the derived 
kinematic parameters will be discussed below. 

For the kinematic fit, the template is first brought to the eff'ec- 
tive instrumental resolution of the observed spectrum by convolv- 
ing with a kernel that is a Gaussian in wavelength. The convolution 
with the LOSVD G is performed in log A (equivalently: velocity), 
and the result is then resampled to the same wavelength grid as 
the data 5, thus avoiding having to resample the data beyond what 
has been done during the data reduction. The functions p''"^\A) and 
q'''^\A) are multiplicative and additive correction polynomials of or- 
der m and n, respectively. These polynomials are needed to correct 
any low-frequency diff'erences between the galaxy and template 
spectra due to insufficient flux calibration, contamination by the 
continuum of the lensed background galaxy and other eff'ects which 
are not related to th e kinematics of the lens galaxy. In contrast to 
ICzoske et al.l (l2008h . where only a short wavelength range around 
a single spectral feature was used and a linear correction function 
was sufficient, we here use the full wavelength range. The coeffi- 
cients are determined by a linear fit nested within the non-linear 
optimization for the kinematic parameters v and 0-^. The orders of 
the polynomials are determined by inspecting large-scale systemat- 
ics in the residuals from the fits; we obtain satisfactory corrections 
for m = 4 and n = 6. 



^ http://www.r-project.org. The slacR package can be obtained 
from the first author on request. 
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Figure 3. Global spectra obtained by adding spectra within elliptical apertures with diameters between 5 and 14 spaxels, matched to the surface brightness of 
the objects. Wavelengths are in A. The three spectra top left are from the pilot programme and cover a different wavelength range then the spectra from the 
large programme. Emission lines from the sources are marked. The narrow residuals of sky lines at 5577 A, 5890 A and 6300 A have been interpolated over 
for display purposes. 



A number of features in the spectra were excluded in the kine- 
matic fits. Possibly present Balmer lines in the lens spectrum are 
due to younger stellar populations that are not adequately described 
by the late- type template stars. The strength of the MgB band and 
Na D lines is enhanced compared to that in the galactic stars used 
as templates (Barth et al. 2002). Also masked were night-sky emis- 
sion lines and atmospheric absorption features, as well as emission 
lines from the lensed background sources. Finally, a 3-cr clipping 
algorithm (three iterations) was applied to detect possible outliers. 

The noise e is assumed to be normally distributed with mean 



zero and wavelength-dependent dispersion cr^. Since vipgi does not 
produce noise estimates that take the details of the reduction pro- 
cedure into account, we have to estimate the noise spectra after the 
fact, using a model including readout noise and photon noise as de- 
scribed in Sect.lUand Appendix lAl The average signal-to-noise ra- 
tio used to characterize a spectrum was determined by dividing the 
spectrum by its noise spectrum and taking the median. Signal-to- 
noise ratios are thus given per wavelength pixel of width 0.644 A. 

The noise was propagated to error estimates on the kinematic 
parameters using a Monte-Carlo method. Gaussian noise was added 
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Figure 4. Kinematic fitting of spectra of system J0935. Top: global spec- 
trum summed in an aperture of radius 2.3 arcsec, S IN ^ AS per wavelength 
bin. Centre: single spaxel with S IN ^ 23 per bin. Bottom: single spaxel 
with S IN ^ 12 per bin. The top panel of each pair shows the spectrum and 
the template fit in red. The bottom panel shows the residuals with the ex- 
pected noise level in green (Appendix |A). The grey areas mark regions that 
were excluded from the fits. Wavelengths are given in the rest frame of the 
lens (z = 0.3475). The [O ii] line from the source can be seen at 4057 A in 
this plots. 



are defined for this purpose as those giving x 



+ 1, where 



X^^^ is the minimum value obtained with the set of stellar spectra. 
The rms values of cr^, obtained with these template candidates lie 
between 5 and 10 % (maximum value is 10.6 % for J0959). For the 
two-dimensional kinematic maps, template mismatch will mostly 
have an effect on the overall level, but not on the structure if stellar 
populations across the galaxy are homogeneous. 



5.2 Results 

Example fits for an aperture-integrated spectrum and spectra from 
individual fibres are shown in Fig. |4] The template fits generally 
reproduce the observed spectra to the expected noise level. 

Fig. |5] shows the maps of velocity, velocity dispersion and 
spectral signal-to-noise. The kinematic maps are restricted to spec- 
tra with an average signal-to-noise ratio S IN > ^ (per pixel) for 
which kinematic parameters can be reliably measured. Two-thirds 
of the sample show little structure in the maps and display kine- 
matics typical for nearly-isothermal pressure- supported slow rota- 
tors. Clear rotation patterns can be discerned in the velocity maps 
for the remaining third of the sample (e.g. J0959, J1251 or J2238). 
Kinematic maps of local galaxies in general show a central peak of 
velocity dispersion; at the redshifts of our sample, the spatial reso- 
lution of the integral-field spectrograph is not sufficient to resolve 
this peak. 

Fig. [6] compares our velocity dispersion measurements on the 
global (aperture-integrated) VIMOS spectra to measurements on 
the spectra from the Sloan Digital Sky Survey (Table[TJ- The agree- 
ment between the different measurements is generally good, which 
gives us confidence in the quality of the spectra and the reliability 
of the analysis method. 

For the SDSS spectra we compare two in dependent measure- 
ments: the values listed bv feolton et al ] 12008) and new values de- 
termined with our code using the same wavelength ranges and tem- 
plate spectra as for the VIMOS spectra. The middle panel of Fig. [6] 
compares these measurements, which differ only in the analysis 
method. The error bars on the slacR measurements include only 
the effect of noise in the spectra and do not take into account sys- 
tematic effects such as choice of template and other parameters. 
The good agreement within the errors shows that systematic effects 
contribute only little to the uncertainties. A notable outlier is J0935 
for which we measure c r - 330kms~^ compared to 396kms~^ 
from lBolton et al.l (l2008h . 

The right panel compares measurements on the SDSS spectra 
to those on the VIMOS spectra, using identical methods. The dif- 
ferent instrumental resolution (R = 1800 for SDSS, R = 2500 for 
VIMOS) has been taken into account. There appears to be a slight 
offset of the velocity dispersions obtained from the SDSS spec- 
tra compared to those from the VIMOS spectra. A similar trend 
is visible in the middle panel, while the compariso n of the VI- 
MOS measurements to the SDSS measurements from lBolton et al.l 
([2008.) does not show any offset. 



to the best-fit model with wavelength-dependent dispersion given 
by the corresponding noise spectrum and the resulting spectrum 
was fitted in the same way as the original spectrum. The quoted er- 
rors on the parameters of the original spectrum are the 16% and 
84% percentiles of the kinematic parameters v and cr^ obtained 
from 300 such realizations. 

The errors due to template mismatch can be estimated from 
the distribution of best-fit values determined with 755 stellar spec- 



6 SECONDARY SOURCES 

Some of the VIMOS/IFU fields in this sample contain objects in 
addition to the primary targets. For these fields we show larger- 
scale cut-outs from SDSS r-band images and label them by number. 
We extracted aperture spectra from the VIMO S data cubes and de- 
termined their redshifts by cross-correlation dTonrv & DavisllT979l) 
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Figure 5. Results of kinematic analysis of the SLACS/IFU sample. For each target, the left plot shows the map of systematic velocity in the lens galaxy and 
the middle plot shows the map of velocity dispersion. All velocities refer to the rest frame of the lens galaxy. The right plots show the mean signal-to-noise 
ration (per pixel of width 0.65 A) of the spectra. The axis labels are in arcsec. The maps are oriented such that north is to the top and east to the left. 
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Figure 6. Comparison of velocity dispersion measurements. The VIMOS spectra used here are summed over all spaxels within t he same aperture as the SDSS 
spectra and analysed using our method (slacR). For the SDSS spectra, values measured with slacR and values obtained from lBolton et alj 1200 8l , BOS) are 
used. 



Table 3. Data for additional sources in VIMOS/IFU fields. Coordinates are 
in arcsec relative to the position of the lens galaxy (always labelled as "1") 
in each field. Redshifts marked with "?" are suggested by peaks in the cross- 
correlation but cannot be confirmed unambiguously by visual inspection. 
Objects with z = are stars. 
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with a set of SDSS template spectra. Coordinates relative to the pri- 
mary target and redshifts are listed in Table [3] The table also lists 
redshifts for the lens galaxies obtained by the same method in order 
to ensure that relative redshifts are accurate. Most of the additional 
objects are consistent with being at the same redshift as the lens 
galaxy. 



7 CONCLUSIONS 

In this paper, we have presented integral-field spectroscopic data 
obtained with VIMOS/IFU on a sample of 17 early-type lens galax- 



ies selected from the SLACS survey. The data permit spatially 
resolved reconstruction of the stellar kinematics in these galax- 
ies, presented as maps of systematic velocity and velocity disper- 
sion. The sample, which spans a redshift range of 0.08 to 0.35 
for the lens galaxies, is wel l- suited for a compar ison with local 
samples, such as SAU RON (lEmsellem et al.ll2004l) and ATLAS 
(ICappellarietal.ll201lh . 

The distinction of slow-r otating and fast-rotati ng galaxies 
found in t he SAURON sampl e (lEmsellem et al.ll2007h and earlier 
work (e.g. IPavies et is also evident from the kinematic 

maps for our sample. About a third of the galaxies in the sample 
show clear evidence for rotation in their velocity maps. 

While the spatial resolution of integral-field spectroscopic 
observations of galaxies at cosmological distances is necessarily 
much coarser than for local galaxies, the present sample oflTers the 
unique advantage that all the galaxies act as gravitational lenses and 
therefore oflfer the possibility to model their mass distribution using 
the two complementary methods of stellar dynamics and gravita- 
tional lensing. A Bayesian implementation of a fully self-consistent 
modelling algorithm for the joint analysis of stellar dyna mics 
and g ravitational lensing was developed by ( Barnabe & Koo pmansI 
|2007|) and applied t o subs ets of this sample bv lCzoske et al.1 (l2008h 
andl fiarnab^ et al .l (l2009l). Modelling a nd analysis of the full data 
set is presented in Barnabe et al.l (l201lh . 
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APPENDIX A: ESTIMATION OF NOISE VECTORS 

Noise vectors can be estimated by following the reduction process 
from the initial raw data through the final data product. Using a 
model that takes into account photon noise (from object and sky) 
and readout noise, it is possible to express the variance in the data 
at each step in the reduction as a simple rescaling of the current 
intermediate data product. 

The raw image frames have values d, expressed in ADU; these 
contain in addition to the useful data the bias b. The noise variance 
is the quadratic sum of the photon noise due to the data values d-b 
and the readout noise r, given in electrons. In terms of electrons, the 
signal (without bias) is {d - b)/g and the variance is thus 

al^=l(d-b) + r^ [e-]. (Al) 

The conversion from electrons to ADU is here given by the gain 
factor g (in ADU/e~). 
Variance scale as 

o-l = c^(t\ (A2) 

when the data are rescaled according to = cf\ . These relations 
preserve the signal-to-noise under rescaling. The variance of the 
raw images in terms of ADU is thus 

o-l^ = g(d-b)+R^ [ADU], (A3) 

where R - gr is the readout noise in ADU. 

The preliminary reduction step in vipgi subtracts the bias 
frame and produces BFC files. Compared to the raw frame, the vari- 
ance is increased by the noise in the combined bias, ^ Since 
the readout noise is small anyway, we neglect this term. 

BFC = d-b (A4) 
cr^p^ = gBFC + R\ (A5) 

The following reduction step applies the wavelength cali- 
bration, rectifies the two-dimensional spectra, and extracts one- 
dimensional spectra from these. The output files are names BFCEx. 



Two-dimensional kinematics ofSLACS lenses -IV. 13 



The EXR2D extension contains the rectified version of BFC. Neglect- 
ing the correlation of pixel noise introduced by the resampling, we 
have, as for the BFC: 



(A6) 



The ID spectra in EXRID are created from EXR2D by the Horne 
optimal extraction algorithm ( Horne 1986; Zanichelli et al. 2005), 
which sums over the spatial direction of the 2D spectrum using an 
estimate of the spatial profile of the light distribution as weights. 
The algorithm permits direct estimation of the noise variance in the 
extracted spectra; this is, however, not implemented in vipgi. We 
replace the Horne variance by the variance of a straight sum over 
Nc = 5 columns in the 2D spectrum: 



EXRID 



^ EXR2D 



= gEXRlD + NcR^ 



(Al) 



(A8) 



Comparison of the noise predicted by Eq. (jASJ to real data in pure 
sky fibres (no object contribution) indicates that the noise is indeed 
somewhat overestimated, although not dramatically so. 

The relative fibre transmissivity ? of a given fibre is estimated 
from the measured flux in a sky emission line and stored in a FITS 
table. Its application results in 



EXRIDTC = 


- EXRID 

t 


(A9) 


2 

<^EXR1DTC ~ 


1 2 

^2 *^EXR1D 


(AlO) 




- - EXRID + ^R^ 

t t f 


(All) 




^ EXRIDTC +-^R^ . 


(A12) 



The flux calibration uses the sensitivity function s{A). It acts like 
any other rescaling. The sensitivity function is normalized to an 
exposure time of 1 s, hence it has to be rescaled to exposure time T 
(in seconds). The estimate of the sky is based on a large number of 
fibers and has negligible eff'ect on the variance. 



EXRIDFLUX = 



1 

g 



(EXRIDTC 
1 



SKYID) 



_ EXRIDTC + 

Tst Ts T^s^r 

Nc 



-^EXRIDFLUX. 



(A13) 
(A14) 
(A15) 



The overline indicates that the sky has not been subtracted. 

The final data cube is generated by combining A^exp exposures 
using the median which results in a somewhat larger variance than 
the arithmetic mean. Simplifying again, we use the variance for the 
latter: 



CUBE 



— y 



EXRIDFLUX 



EXRIDFLUX 



(A16) 



(A17) 



Note that due to the dithering between exposures some of the spax- 
els in the cube arise from a combination of fibres from diff'erent 
quadrants which have diff'erent gain g and sensitivity functions s. 
Assuming that we can use some eff'ective values averaged over the 
four quadrants, the factor gl st'm. the second row of Equation ( IA17I ) 
can be taken out of the sum. This way, the variance of the cube 
can be obtained from a mock cube CUBE which is produced in the 
same manner as the actual cube CUBE, only with sky subtraction 
turned off". 

All the parameters in Eq. dAlVl l are stored as FITS header key- 
words or in table extensions of vipgi products. 

In the final step, the global spectrum is created by summing 
over A^fib spaxels (defined by an aperture) in the final data cube, 
giving 



SPEC 



= cube 

~ Yj ^cube 
N..„ Tst ^ 



A^exp Tst 



SPEC + 



(A20) 



(A21) 



(A22) 



(A23) 



Equations (IA20I) and (IA21b provide the recipes for the final 
spectrum. The readout noise r is written in terms of electrons, as it 
usually is given in the fits headers. 



